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Abstract: We introduce a novel technique designed to look for signatures of new physics 
in vector boson fusion processes at the TeV scale. This functions by measuring the polar- 
ization of the vector bosons to determine the relative longitudinal to transverse production. 
In studying this ratio we can directly probe the high energy £' 2 -growth of longitudinal vec- 
tor boson scattering amplitudes characteristic of models with non-Standard Model (SM) 
interactions. We will focus on studying models parameterized by an effective Lagrangian 
that include a light Higgs with non-SM couplings arising from TeV scale new physics as- 
sociated with the electroweak symmetry breaking, although our technique can be used in 
more general scenarios. We will show that this technique is stable against the large uncer- 
tainties that can result from variations in the factorization scale, improving upon previous 
studies that measure cross section alone. 
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1. Vector Boson Fusion as a Probe of New Physics 

The Large Hadron Collider (LHC) was built to elucidate the physics behind electroweak 
symmetry breaking (EWSB). In a sense, it must succeed in finding some new physics be- 
cause the partial wave amplitudes for VlVl — ► ViYl scattering, 1 calculated in the absence 
of a Higgs or other new physics, begin to violate unitarity at the TeV scale. Therefore, 
either new weakly-coupled light particles must come in to unitarize the amplitudes, or we 
will see new strong interactions in the electroweak sector. 

While many models of EWSB have been proposed, precision experiments such as LEP 
seem to favor a model employing a 0(100) GeV scalar with the quantum numbers and 
approximate couplings of the Standard Model (SM) Higgs [1, 2]. Many models of new 
physics already include such a particle, oftentimes with couplings deviating slightly from 
those of the SM, e.g. little Higgs [3] and holographic Higgs models [4]. Ideally, such models 

1 By Vl we denote a longitudinally polarized electroweak vector boson. 
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would be identified and studied at the LHC through the production of their intrinsic new 
particles. However, the finite energy reach and large backgrounds at the LHC could make 
discovering any new states very difficult. 

Thus we will focus on these non-SM light Higgs scenarios, both because they are 
favored by precision data and because they are perhaps the most difficult to distinguish 
from the SM. To study these setups we will take a model-independent approach, employing 
an effective field theory to parameterize the effects of new physics [5, 6, 7, 8, 9, 10]. We 
will see that the general phenomenology of the Higgs sector is captured by the coefficients 
of a small number of dimension-6 operators [11, 12], only one of which is relevant to the 
vector boson fusion process we wish to study. 




Figure 1: Illustration for vector boson fusion. 



Vector boson fusion (VBF) is the process in which vector bosons radiated by initial 
state quarks scatter into vector bosons (see Fig. 1). This process is intimately tied to 
EWSB: just as the pion is a Nambu-Goldstone boson (NGB) and irn scattering can be used 
to understand chiral symmetry breaking, at high energies longitudinally polarized vector 
bosons take on the behavior of the NGBs from EWSB. In the absence of a Higgs boson 
or other new physics responsible for the EWSB, the scattering amplitudes probed by VBF 
would violate perturbative unitarity [13, 14, 15, 16] at around 1 TeV (see the discussion in 
appendix A). Furthermore, if the Higgs boson does not have the exact couplings to vector 
bosons as predicted by the SM, then the necessary cancelations will not occur and one will 
still observe an E 2 growth in the amplitudes until new physics comes into play. It is by 
measuring this growth that we can hope to observe the effects of physics beyond the SM, 
even in scenarios where we only see a light Higgs-like particle [12, 17]. 

In this article we will introduce a novel technique designed to analyze VBF processes 
and observe the E 2 growth in longitudinal gauge-boson scattering amplitudes mentioned 
above. We will begin by introducing our notations and framework in Section 2. To mo- 
tivate our new technique, in Section 3 we will update past analyses of VBF (specifically 
[18] and [19]) taking into account the effects of parton showering and jet clustering. We 
will show that these analyses, which infer the E 2 amplitude growth from cross section in- 
creases, carry large 0(100%) uncertainties due to factorization-scale ambiguities that affect 
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jet tagging and vetoing. 2 Then, in Section 4 we will describe our technique designed to re- 
duce these uncertainties by measuring the relative production of transverse to longitudinal 
modes, focusing on the fully reconstructable semi-leptonic decay of the VV system. We 
will demonstrate that this measurement is sensitive to anomalous Higgs-gauge couplings 
while at the same time being robust against the scale uncertainties that challenge cross 
section measurements. In Section 5 we will discuss potential improvements in the analysis 
and comment on other applications of the technique. Section 6 contains our conclusions. 

2. Theoretical Setup 

In the formulation of a general effective theory of the SM-like Higgs sector [5, 6, 11] most 
of the operators are tightly constrained [7, 8, 9, 10] because of their otherwise excessive 
contributions to the electroweak observables, such as the /^-parameter, oblique parameters, 
and triple gauge boson self-interactions. There are only two dimension-six operators that 
are genuine interactions in the Higgs sector not subject to the current experimental con- 
straints, d^{H^H)d^{H^H) and (H^H) 3 . We note that as both operators are composed 
from the singlet operator H^H they may serve to probe not only EWSB physics, but also 
other physics beyond the SM. For a given theoretical framework, the coefficients of these 
operators may be calculable [12], and by measuring them we can hope to learn about any 
new physics. Even in some strongly coupled models for which these may not be calculable, 
the measurement of a non-zero value can give important clues to the structure of new 
physics. Now, the second operator above does not have derivative couplings, so its effect 
on the behavior of the Vl scattering amplitudes at high energies should be sub- leading [11]. 
We therefore focus on the former and parameterize it (following [12]) as 

C D ^{rfH)d^H), (2.1) 

where the coefficient ch is naturally of 0(1 — 4ir) depending on whether the underlying 
theory is weakly or strongly coupled, and / is the characteristic scale of new physics, 
typically expected to be round Attv if the new physics is associated with EWSB. 

Upon expanding around the electroweak VEV v, this operator contributes terms which 
add to the kinetic terms of H. After imposing canonical normalization on the fields, 
the result is a modification to the Higgs couplings. Ref. [12] parameterizes the resulting 
modified Higgs-gauge coupling in the zero-momentum limit as 

where £ = v 2 /f 2 . This modified coupling has important phenomenological consequences 
because it deviates from the SM prediction. At high energies and for |c#£| > 0.1 3 this 

2 While one might expect to eventually overcome these uncertainties through a combination of data- 
driven calibration and higher order calculations, the analysis we describe herein avoids them altogether. 

3 For smaller values of \ch£\ the dominant non-SM effects enter as interference terms proportional to csr£ 
rather than (ch£,) 2 - Also, in this case the anomalous energy dependence of the longitudinal cross section 
goes as E 2 instead of E 4 . 
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modification leads to an incomplete cancelation in the amplitude for longitudinal vector 
boson scattering and the cross section grows as 



a{V L V L -» V L V L ) « (^e) o{VlVl -> VL^)no-higg S - (2.3) 

which can be seen by considering the NGB scattering as shown in appendix A. In what fol- 
lows we will study means of measuring this behavior. Note that, as discussed in appendix A, 
the W^W£ scattering amplitudes calculated in this framework violate perturbative uni- 
tarity when 

This is the point at which we expect new physics to come into play. In what follows we 
will limit our analyses to 

^fs^v < 2 TeV. (2.5) 

This corresponds to a coupling value \ch£\ ~ 0.6. We will take this as an upper limit for 
our analyses. Of course, looking beyond this energy range would be interesting and should 
be attempted at the LHC, but any deviation from the SM expectation would no longer 
carry the same effective Lagrangian interpretation. Also, note that for larger couplings and 
lower scales of new physics some higher dimensional operators could become relevant and 
it would be more appropriate to think of the used in our analysis as parameterizing a 
new physics form factor, rather than as the coefficient of a particular operator. 



3. Leptonic and Semi-Leptonic Channels Revisited: Scale Uncertainties 

The most straightforward way to probe the behavior of Eq. (2.3) would be to measure the 
resulting increase in VBF cross section at higher energies. This is a well studied topic, 
with many different analyses having been performed (see, for instance, [18, 19, 20, 21, 22, 
23, 24, 25, 26, 27, 28, 29, 30, 31]). Here we will revisit these analyses taking into account 
the effects of parton showering. We will see that while the cuts developed in past analyses 
remain essential in suppressing background and isolating VBF signals, one encounters large 
factorization scale uncertainties that affect rate measurements and must be overcome to 
detect new physics in VBF. 

Sophisticated acceptance cuts have been developed to suppress the SM background 
and isolate the longitudinal gauge boson scattering in VBF processes at high energies. It 
has been a common practice to impose a high px cut on the reconstructed gauge bosons 
or their decay products, require one or two forward (backward) energetic jets, and demand 
that the central detector region remain relatively free of hadronic activity. The first few 
cuts ensure that we observe hard scattering processes with the gauge bosons emitted by 
energetic quarks [32, 33, 34], while the last cut is designed to reduce background by taking 
advantage of the fact that VBF is a purely electroweak process with no color exchange [35] 
and VlVl scattering tends to produce fewer central jets than other electroweak processes 
[21]. Using this sort of cut, it was concluded [18, 19] that reasonable sensitivity can 
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be achieved for TeV scale strongly interacting new physics at the 14 TeV LHC with an 
integrated luminosity of 100 fb _1 . 

We revisit the WW analyses with the theoretical framework as discussed in the previ- 
ous section. We consider both fully leptonic [18] and semi-leptonic [19] decays of the vector 
bosons. For the sake of illustration, we concentrate on the W + W~ final state. Our VBF 
parton-level results are generated using the full 2 — ► 6 matrix element for 



without making the effective W approximation [36, 37, 38]. In so doing, wherever appro- 
priate, we have included other 0(oft EW ) processes as background to the channels. Our 
PDFs are those of MRST2004 [39]. To generate the jet-level samples we shower parton- 
level results using Pythia 6.4.21 [40] with a virtuality ordered shower, cluster the visible 
final state particles into 0.1 x 0.1 y — 4> cells between — 5 < y < 5, and produce R = 0.7 
anti-fcr [41] jets using Fast Jet [42]. To sample PDFs and shower our results we must choose 
a factorization scale for the gauge boson scattering processes. The natural choice of the 
factorization scale is of the order of m w , with corrections from the pt of the scattering 
quarks. We parameterize the choice of scales via 



where (5 is an 0(1) parameter. 

We begin by adopting the selection cuts of [18] to study the fully leptonic W + W~ final 
state, as detailed in Table 1. Using these cuts, we calculate the parton level cross sections 
for a light Higgs scenario 4 with various anomalous couplings parameterized by c H £. The 
parton level results for a few representative scale choices of (3 are listed on the left-hand 
side of Table 2. They are consistent with those of [18]. At this level, the uncertainty in 
rate is only around O(10%), which can be attributed entirely to the PDFs. When we 
include showering, hadronization and jet clustering, the scale \x dictates the virtuality at 
which the parton shower is started, in addition to controlling the sampling of PDFs. On 
the right-hand side in Table 2, we present the cross sections for the showered and clustered 
results with a few representative scale choices. We see that the uncertainties can now be 
as much as O(100%). This is because small changes in p, result in large changes in the 
behavior of the associated forward jets. A higher value of /i could lead to harder radiation 
that will sink forward jets below the tagging criteria, or it could lead to the parton-shower 
emission of a veto jet. As the uncertainties from varying the scale {(3 = 0.5 — 2.0) would 
normally set the systematic theoretical errors, such large uncertainties in rate would make 
it difficult to distinguish the presence of anomalous couplings, even for large values of cj/£. 

We next explore the situation for the semi-leptonic mode of W + W~ decay. We employ 
cuts inspired by [19] as shown in Table 3. The results of this analysis are shown in Table 4, 

4 Here and henceforth, we will take a light Higgs boson mass as mu = 100 GeV for illustration and for 
comparing with the early studies in the literature. This will make no numerical difference with other rriH 
values as long as it is well below the 2m w threshold. 



qq' -» qq'W + W 



qq' //', 



(3.1) 




(3.2) 
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Leptonic Cuts 


Jet Cuts 


\y(l)\ <2.0 
pr(l) > 100 GeV 
A PT (ll) > 440 GeV 

cos 4>u < —0.8 
M(ll) > 250 GeV 


^(j'tag) > 0.8 TeV 
3.0 < |j tag | < 5.0 
PrO'tag) > 40 GeV 
Mivcto) > 30 GeV 
|y(jVcto)| <3.0 



Table 1: The cuts of [18] for the leptonically decaying W + W final state. The signal selection 
requires that we tag at least one jet (jtag) and to veto extra central jets (jvcto)- 





Parton Level [fb] 


Jet Level [fb] 




P = 0.5 


p = 1.0 


= 2.0 


P = 0.5 


= 1.0 


= 2.0 


0.4 


0.015 


0.013 


0.012 


0.015 


0.009 


0.005 


0.2 


0.013 


0.011 


0.010 


0.013 


0.006 


0.004 


0.0 


0.011 


0.090 


0.008 


0.012 


0.007 


0.004 



Table 2: Cross sections [fb] for VBF with W + W~ final states decaying into e and fi for various 
anomalous Higgs-gauge couplings and at different factorization scales parameterized by Eq. (3.2). 
The cuts used to generate these results are those of [18]. The set of cross sections on the left are 
computed at parton level, while those on the right correspond to results after the parton shower 
and hadronization. 



again demonstrating a relatively stable signal at the parton level (left-hand panels) and an 
O(100%) uncertainty at the jet level (right-hand panels). As with the fully leptonic system 
considered above, the large uncertainty is once again attributable to the parton-shower 
treatment of the forward jets using different scales. 

We find that VBF cross sections with complicated kinematical cuts are extremely 
sensitive to the exact scale choice one uses, and it is insufficient to simply characterize the 
hard scattering with a rough estimate of fj,. While we have only considered the effects of 
varying \x on signal rates, the background is also susceptible to these uncertainties. Even 
without considering the effects of the veto cut, which can only increase sensitivity to n, it's 
natural to associate 0(2 x) K- factors with high multiplicity QCD events. Therefore, before 
drawing any conclusion about the presence of new physics one would have to understand 
these systematics. In principle, the theoretical uncertainty may be reduced through higher 
order calculations that can give us a better idea of the appropriate scale choice. Substantial 
efforts would be needed both in theory and in experiments before to bring this uncertainty 
under control. With this in mind, in the next section we will present a new tool to 
circumvent the difficult issue of the factorization-scale dependence. 
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Pass conditions 


Veto conditions 


£(jtag) > 300 GeV 
2 < |y(jta g )| < 5 
PTO'tag) > 20 GeV 

PT (^ recon .) > 320 GeV 

|y(W had )| <4 


PT(j'mini) > 25 GeV 
|y(j'mini)| < 2 

130 GeV < mvt/j < 240 GeV 



Table 3: W + W~ semi-leptonic decay cuts inspired by [19]. These require two tagged jets (jtag) 
and two reconstructed Ws (Wrocon.)- If the events contain two soft, central jets (jmini) they are 
vetoed. The cut on the jet-TV invariant mass is designed to reduce top quark backgrounds. 





Parton Level [fb] 


Jet Level [fb] 




P = 0.5 


= 1.0 


[3 = 2.0 


= 0.5 


= 1.0 


[3 = 2.0 


0.4 


0.95 


0.81 


0.73 


0.53 


0.38 


0.26 


0.2 


0.82 


0.72 


0.64 


0.43 


0.33 


0.24 


0.0 


0.73 


0.64 


0.57 


0.40 


0.29 


0.21 



Table 4: Cross sections [fb] for VBF with the semi-leptonic decay of W + W~ for various anomalous 
Higgs-gauge couplings and at different factorization scales parameterized by Eq. (3.2). The cuts 
used to generate these results are those of Table 3. 

4. Polarization Measurements 

With the uncertainties detailed above as our motivation, we propose a new technique to 
probe the anomalous couplings in a robust way. Our basic idea is to look for the rela- 
tive increase in longitudinal vector boson production by comparing it to the production 
of transverse modes. Unlike the overall cross section, which is sensitive to the behavior 
of the forward jets, the relative transverse to longitudinal production rates should be sta- 
ble against different scale choices because it depends only on the VV — > VV scattering 
amplitude. 

To measure the polarization of a vector boson we need to reconstruct the four-momenta 
of its decay products and measure their distribution with respect to a polarization axis. If 
one chooses the polarization axis to be the gauge boson direction of motion (Fig. 2), then 
a simple spin-analysis predicts that in the V rest frame the transverse and longitudinal 
polarizations will be distributed as 

P±(cos0*) = 3 (l±cosr) 2 , P L (cos6*) = |(1 -cos 2 ^) (4.1) 
8 4 

where 6* denotes the angle between the parton and the gauge boson direction of motion 
in the gauge boson rest frame. 
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Figure 3: The distributions of cos 8* for different anomalous couplings at parton level (left) and 
for fully showered, hadronized, clustered, and reconstructed events (right). All distributions are 
normalized to the same area. 



V 




I 



To measure these distributions experimen- 
tally, we need to fully reconstruct the gauge 
boson pair center of mass and each gauge bo- 
son's direction of motion in this frame. To ac- 
complish this we will focus on the semi-leptonic 
decay channel of the V V system as this allows 
full reconstruction of the system while mini- 
mizing the SM background by requiring lep- 
tons and missing energy. The semi-leptonic 
channel also significantly increases the signal 
event rate. For this we will rely upon jet sub- 
structure techniques to reconstruct the hadron- 
ically decaying gauge boson [19]. We will focus 
on studying the W + W~ final state, although 
we will take into account the background from 
other VBF processes like W^VF 1 * 1 and W^Z that enter because we can not distinguish 
the sign of a hadronically decaying vector, nor can we always distinguish a hadronically 
decaying W from a Z. Later in this section we will comment on the SM 0{a 2 s ) and 0{a\) 
backgrounds. 



Figure 2: The polarization axis used to 
measure 6*. Note that this is measured in 
the rest frame of the W, and the W direction 
of motion is defined with respect to the WW 
center of mass. 



4.1 Leptonic Polarization 

We begin with the polarization analysis for the leptonic side of the decay. We first study 
the parton-level results, then we will turn on the full simulation (parton-showering and jet 
clustering) to see that they are largely unchanged. The cuts we use are those in Table 3. 
Before proceeding further, we encounter a subtlety in the reconstruction of the leptonic 
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Leptonic W 


Hadronic W 






f p 

JL 


Jl 


f P 
Jl 


Jl 


a [tbj 


-u.u 


0.77 


0.74 


0.71 


0.55 


o.oo 


-0.4 


0.58 


0.52 


0.49 


0.40 


1.12 


-0.2 


0.33 


0.30 


0.23 


0.24 


0.60 


0.0 


0.27 


0.25 


0.17 


0.22 


0.62 


0.2 


0.34 


0.31 


0.24 


0.26 


0.65 


0.4 


0.42 


0.39 


0.32 


0.32 


0.73 


0.6 


0.46 


0.42 


0.40 


0.38 


0.87 



Table 5: The fraction of longitudinally polarized vector bosons for different anomalous couplings 
at parton level and jet level //, reconstructed in hadronic and leptonic decays. Also listed are 
the jet-level cross sections. These results are after the cuts of Table 3. 

system: While the neutrino four-momentum is constrained by the on-shell W condition, 
it is only determined up to a discrete ambiguity. One finds two candidate four-momenta 
at the same azimuthal angle but separated from the charged lepton by a fixed rapidity 
difference. In what follows we will simply use the average cos 9* value from both solutions 
as an approximation of the true value. This is acceptable because we are working in 
a boosted regime where the difference in rapidity between neutrino and lepton is small, 
making the curvature effects from the (y, (f>) system sub-leading. The resultant distributions 
are shown in Fig. 3, at parton level (left panel) and after the hadronization (right panel). 
The characteristic shapes with different couplings are quite distinctive. In Table 5 we 
compute the cross section for each anomalous coupling and fit it to the transverse and 
longitudinal distributions of Eq. (4.1) using 

P(cos0*) = /zPl(cos0*) + /+P + (cos0*) + /_P_(cos0*) (4.2) 

where the P are normalized probability distributions of cos 6* and the / are subject to the 
constraint ^ / = 1. As one can see from comparing the jet and parton level figures, the 
results are remarkably stable under a full simulation. 

In Figure 4 we plot the projected event distributions and associated statistical errors 
both for the SM and for an anomalous scenario with ch£, = —0.4, given 100 fb" 1 of 
luminosity. The shape difference between the two samples is clearly visible. To estimate 
the luminosity necessary to probe a given coupling, one can use that the signal scales 
roughly as (c#£) 2 , as discussed before. However, the precise reach of the LHC in discerning 
anomalous couplings will require a more thorough accounting of background (we discuss 
this further in Section 4.4). Further, we have not made an effort to optimize the statistical 
power of the analysis and there are other channels that contribute to the signal, such as 
W + W + ,W^Z and ZZ. In addition, one can extract more information from each event, 
as we will now see. 
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Figure 4: Projected distribution and associated statistical uncertainties of cos 9* for the leptoni- 
cally decaying vector using 100 fb _1 of luminosity. 

4.2 Hadronic Polarization 

It is possible to further improve the discriminating power of polarization by considering 
both sides of the VV system together; by looking for the expected correlation between 
both states one can hope to gain additional discriminating power. 

To see the correlation effect, consider Fig. 5, which shows the parton-level cos 9* dis- 
tributions for both sides of the VV system in SM and non-SM scenarios. For now, we plot 
cos 9* on the hadronic side for the down- type quarks. In the non-SM scenario we see a 
rapid rise in the central region of the plot near cos 9* ~ 0. This indicates that the results 
are correlated; when we see a Vl it is likely to be accompanied by a Vl because only the 
VlVl final state sees the E 2 growth characteristic of with non-SM effects. In practice the 
situation is slightly more complicated because we cannot label the light quark states once 
they shower and hadronize (e.g. we cannot distinguish a u from ad), so the distributions 
we measure are symmetrized. However, the distributions still carry additional discriminat- 
ing power, as one can see from the distributions in Fig. 6 and Fig. 7, and Table 5. Note 
that in fitting the symmetrized distributions we only fit to data from < | cos#*| < 0.7. 
In the regime where | cos 9*\ > 0.7 one subject becomes very soft and the technique breaks 
down (although, of course, the leptonic analysis still works here). 

To perform this analysis we had to look at the hadronically decaying V using subjet 
techniques (for a short overview of jet algorithms and their behavior, see appendix B). In 
particular, we used the kx algorithm [43, 44] with R = 0.25 to cluster the constituents of 
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Figure 5: Distributions of cos (9* at parton-level for both sides of the VV system (labeled with 
subscript H and L for hadronic and lcptonic decays, respectively). The plot on the left is for the 
Higgs with SM couplings, while the one on the right is for cjj£ = —0.6. The scale is individually 
normalized for each plot, going from violet to red as the concentration of events increases. The 
scaling of the color gradient on the right side of each plot is linear. 

each hadronically decaying gauge boson, using the two most energetic subjets (as measured 
in the VV center of mass frame) for our analysis. Note that rather than identifying our 
subjets through a C/A [45, 46] or fey-like unwinding [47, 48], we used fixed small cones (i.e. 
small R). Otherwise, the subjets encompass a large area and become more susceptable 
to contamination from initial state radiation, multiple interactions, and event pileup. The 
choice of a small cone seems to result in a better reconstruction of events, especially at 
high values of cos 9* when there is a large difference in the subjet pts. Furthermore, we 
use fey rather than anti-fey to form our subjets because it more accurately reconstructs the 
softer jet in situations where the jets are nearly collinear (see appendix B). 

One important thing to consider in the subjet analysis is that the results are not 
as robust in going from matrix-element to parton shower as were the leptonic results; 
the curves change shape (compare the parton and jet level results for both sides of the 
decay in Table 5). This is because the diffuse nature of the subjets makes them difficult 
to resolve when they become collinear and/or soft. We note, however, that at the LHC 
we can expect to calibrate subjet measurements for boosted hadronic Ws with large SM 
samples, and while the parton-level to jet-level results may vary, the correspondence should 
eventually be well understood. Thus the leptonic gauge boson analysis is likely to be the 
first tool used, but the hadronic analysis can be added later on. 

4.3 Scale Uncertainty 

We now explore the sensitivity of our results to the choice of the factorization scales. 
As before, we vary the scale choice as parameterized by Eq. (3.2). We present our fully 
showered results for the leptonically decaying W in Table 6 using different scale choices. 
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I cos 9 I I cos I 

Figure 6: Distribution of | cos 8* | at different anomalous couplings for hadronically decaying Ws 
using parton level samples (left) and fully showered, hadronized samples (right). Note that the 
distributions differ more at high values of cos 8* because this is the region in which one jet is 
relatively soft. 





Longitudinal Fraction 




= 0.5 


(3 = 1.0 


(3 = 2.0 


0.0 


0.25 


0.26 


0.25 


0.2 


0.33 


0.33 


0.33 


0.4 


0.40 


0.40 


0.41 



Table 6: The reconstructed fraction of longitudinally polarized leptonically decaying vector bosons 
measured for various anomalous couplings and at different factorization scales. 

One can see from these results that the relative longitudinal to transverse fraction is a 
robust measure of anomalous couplings, steady across different values of (3. The reason for 
this is clear: The showering and hadronization processes of the forward jets are independent 
of the behavior of the polarizations and decays of the gauge bosons. Whatever happens 
with the forward jet-tagging and the central jet-vetoing is isolated from the polarization 
measurements of the final decay products, hence the robust behavior. 

4.4 Background 

In our matrix element calculations, we have included the irreducible electroweak VVjj 
backgrounds. A full treatment of the other background from QCD and top quarks is 
beyond the scope of this work, but here we give some rough estimates and discuss how 
background processes might have an effect upon the calculations. With the sophisticated 
acceptance cuts as adopted here, estimates of the background generally place it at or below 
the signal level. For the semi-leptonic modes of our interest in [19], whose cuts we mimic, 
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Figure 7: Jet-level distribution of cos 9* (labeled with subscript H and L for hadronic and leptonic 
decays, respectively) for the SM Higgs (left) and Higgs with cjj£ = —0.6 (right). The scale is indi- 
vidually normalized for each plot, going from violet to red as the concentration of events increases. 
The scaling of the color gradient on the right side of each plot is linear. 

the background is calculated to be roughly one quarter of the signal rate, so we expect 
the enhanced signal distributions to be discernible. We would like to emphasize that the 
robust feature of the polarization measurements should still persist when comparing with 
the background events, in which the leptonic decaying W is mostly transversely polarized, 
and the di-jets from QCD have no particular angular preference. Thus further separation 
of the backgrounds may be possible by extending this technique to the two-dimensional 
angular plots where one looks for correlations as shown in the last section. 

5. Future Directions 

Here we discuss some further improvements one could make in reconstructing the hadron- 
ically decaying gauge bosons, along with some additional applications and extensions of 
our technique. 

5.1 Improving Hadronic Reconstruction 

Comparing the columns of Table 5, one can see that the polarization of the hadronicaliy 
decaying gauge boson was reconstructed less accurately than the gauge boson that decayed 
leptonically. This is partly because at higher boosts the gauge boson decay products 
become more collinear and are consequently more difficult to resolve individually. While 
it is tempting to simply exclude this region from the analysis, the highly boosted region is 
precisely where we expect to see the greatest effects of new physics, so it is worthwhile to 
try to improve reconstruction in this regime. 

One approach for better resolving the subjets would involve letting the subjet radius 
scale as l/pr, so that the size of the subjets is adjusted to naturally account for their 
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kinematics. This sort of scaling was explored by [49], although in a non-subjet context. 
We note that it is likely one would have to adapt the 1/pr behavior to account for the finite 
calorimeter segmentation and crowded subjet environment if this approach were adopted. 

Another approach, along a different line, would involve abandoning subjet construction 
at very high boosts in favor of an energy splitting technique. Such a technique would split 
the jet from a hadronically decaying gauge boson into two pieces and consider the degree 
of energy sharing between them using something akin to the z- variable of [48] : 

z = mm(E A ,E B )/E tot . (5.1) 

For cos 8* ~ 0, where both daughters are transverse to the V direction of motion we would 
expect z ~ 0.5, while for cos#* ~ 1 we expect z ~ 0. We avoided this technique earlier in 
our polarization analysis because it is not boost invariant, and an alternative definition of 
z using pt instead of E would still be sensitive to the orientation of the V decay. However, 
these effects go as my/ /pr, and so are marginalized as we consider higher boosts. Therefore, 
it might become advantageous to transition to a z variable in the highest px ranges. 

5.2 Additional Applications 

Throughout this article we have focused on scenarios with a light, SM-like Higgs because 
they can be very difficult to distinguish from the SM, making VBF analyses particularly 
useful. However, the VBF polarization measurements we have discussed could be of use in 
studying other models of EWSB. 

Consider, for example, the class of models where VlVl scatter amplitudes are unita- 
rized by the exchange of new spin-1 states (see Fig. 8). Models in this class include the 
so-called Higgsless models, where higher dimensional gauge bosons unitarize the longitudi- 
nal amplitudes, and technicolor, which employs the technirho for this purpose [50, 51, 52, 
53, 54, 55, 56]. 

Phenomenological studies of these models [57, 58] often study the VBF production of 
the new resonances because the heavy spin-1 states involved in restoring unitarity must 
couple strongly to the longitudinal modes of the electroweak gauge bosons. While it would 
be striking to see a bump in the VV invariant mass to reconstruct a resonance, seeing 
this correlated with a measured increase in the longitudinal production would give a true 
smoking gun signal. Furthermore, if the new states have a large width and can not easily 
be seen as a resonance, their presence might be inferred by a broad increase in longitudinal 
mode production near the particle mass. This would also be useful in probing a scenario 
with a heavy Higgs. 

Furthermore, in models of EWSB where a heavy Z' mixes with the SM Z, there are 
even more dramatic effects one could look for. In these models the VBF aptitudes grow 
as E 4 up until the scale of the Z' [59] . These models can be difficult to probe if the Z' is 
too heavy to be produced, but one could in principle observe the unique amplitude growth 
characteristic of these scenarios. 

Another extension of this sort of analysis using VBF as a robust probe of new physics 
would involve making use of azimuthal angle correlations. This sort of analysis has been 
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performed with an eye toward distinguishing the spin of the unitarizing particle [60]. It 
would be interesting to see if this could yield an additional handle on background suppres- 
sion or signal enhancement. 

Finally, we note that while our article focuses on the E 2 growth in VlVl amplitudes 
associated with new physics, for many interesting cases the rates are so low that with the 
LHC we can probably only hope to observe this growth when integrated over all higher 
energies. A luminosity upgrade to the LHC [61] might allow for the differential detection 
of this growth: given enough data, one could bin longitudinal fractions according to the 
s of the reconstructed VV system, allowing for a differential detection of the amplitude 
growth. 



by the current electroweak precision data. However, this type of model is the most difficult 
to distinguish from the SM, especially if the new physics particles are very heavy. However, 
it is also the scenario in which VBF can be most useful, because for such a scenario the 
amplitude for VlVl - > VlVl scattering exhibits a non-SM E 2 growth until new physics 
comes into play. 

Past analyses designed to measure this E 2 growth were reviewed and updated to 
account for the effects of the parton shower and jet clustering. While the cuts pioneered 
by these works can be very powerful in reducing the SM background, we demonstrate 
that there is still a significant 0(100%) rate uncertainty attributable to factorization scale 
ambiguities. Thus, we show that in the absence of higher order calculations that might 
give us some guidance on the correct scale treatment, rate information alone may not be 
sufficient to distinguish the signs of new physics. 

We then propose our new technique, which uses the semi-leptonic decay mode of the 
VV system to fully reconstruct events and obtain the decay angle distributions for the 
V daughters. These distributions can be decomposed into longitudinal and transverse 
components, allowing us to measure the E 2 growth in scattering amplitudes associated with 
new physics by looking for the relative increase in longitudinal production. We demonstrate 
that these results are insensitive to the scale ambiguities that trouble rate measurements. 

In closing, we wish to reiterate that polarization measurements of VBF final states 
are a powerful, robust probe of new physics associated with the EWSB. Although we have 
only employed them here to study light SM-like Higgs scenarios, they would be useful in 
more general scenarios of EWSB as long as the longitudinal gauge bosons are significantly 




6. Conclusions 



In this article we have introduced a powerful new 
technique for identifying signs of beyond the SM 
physics associated with the EWSB by probing VBF 
processes at the LHC. 



Figure 8: Unitarity restoration through 
a spin-1 particle. 



We began by motivating our decision to study 
models of EWSB employing a light Higg-like par- 
ticle with couplings deviating from those of the 
SM. Theories with a light Higgs boson are favored 
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involved. Such measurements may prove to be our best tool in understanding the physics 
of EWSB at the LHC. 
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A. Scattering Amplitudes for Longitudinal Gauge Bosons and Partial 
Wave Unitarity 

For completeness, we will here review the high energy behavior of longitudinal gauge boson 
scattering. This will demonstrate why we expect the increase in the scattering amplitudes 
for non-SM Higgs couplings. It will also help us establish the partial wave unitarity bound 
for longitudinal gauge boson scattering, which is of practical importance for our simulation. 
We will make use of the Goldstone equivalence theorem, which says that the scattering 
behavior of the longitudinal gauge bosons is the same as that of the eaten Goldstones, up 
to corrections of order 0(m\y/E). Note that while we will only explicitly calculate the 
behavior of W^W^ — ► W£W£, the other longitudinal gauge boson scattering processes 
are similar. 

We begin with the Lagrangian for the SM Higgs doublet with the additional dimension- 
6 operator we wish to study: 

£ = 1 Tr (m-Cdn + ii 2 rCrdj - ^ Tr (rOrdj 2 + 9 Tr (w f w) d Tr (tCh) (a.i) 

where H = ■ & for real fields 0j (i = 3) and a = (l,cr). The SM Higgs potential 
corresponds to c H = 0. Expanding around the minima (0o) = p/VX one finds new deriva- 
tive interactions proportional to ch- Those relevant to 0+0- — > 0+0- scattering at lowest 
order are: 

C D -v\hcp +( j)_ - ^(f> 2 + 4?- + {f+{d<t>-) 2 + 0+0-50+50- + 2^0+,90-d/i) +h.c. (A.2) 

where we have denoted the shifted 0o field by h and written 0^2 in terms of their charge 
eigenstates 0+. Also, note that there is an additional kinetic term for h: 

£ D f^) 2 (A.3) 

so that in going to canonical normalization we must insert a factor of N = 1/ y/1 + c H v 2 /f 2 
for every h encountered at a vertex. The tree level amplitude becomes 

2 



.ch iN 2 ( c H v 



A^(0+0_ -> 0+0-) = -MX + A - - ( 2Xv + -^s ) +{s^t) (A.4) 
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where rrih = y2\v. Working in the limit s,t^> m h we find 



M{4>+<!>- -» <M_) « (l - ^ 2 ^f j ( 8 + t) (A.5) 

which shows the -E 2 growth in the amplitude that we expect. In this limit, the J = 
partial wave is: 

«=^jC?( i -<£)<- + <>*=i^( i -<£) (a - 6) 

Partial wave unitarity is violated when |Re(ctj)| > 1/2, so the unitarity bound is saturated 
when 

16irv 2 

W " chZ(1-chZN*/4) {AJ) 
To stay clear of this limit, we limit ourselves to studying events for which s < 2 TeV 
(corresponding to \ch£\ < 0.6). 

We note that one may gain further intuition into the longitudinal gauge boson system 
by considering the parameterization 

U = (v + h) e**- B l v (A.8) 

Here we have shifted our field definitions so that the tt transform non-linearly. In this 
language, the relevant terms in the Lagrangian become 



C D —hd-n+d-K- + -\ (vr 2 (d7r_) 2 - ^ + 7r_9vr + avr_) + h.c. (A.9) 

As before, the kinetic term of h is shifted, so we must add a factor of N at every point 
we encounter an h at a vertex. Note, however, that in this case all of the operators come 
with dir terms. Computed in this way, the amplitude 7W(7r + 7r_ — ► 7r + 7r_) shows the same 
behavior as Eq. (A. 4), as it must, but this is the result of a non- cancelation of derivatives 
between the four-point operator and the /i-exchange in the t & s-channels, rather than 
because of a new vertex. 

Using these results we can compare the scattering in the Higgsless case to that of 
the case where the Higgs has anomalous couplings. To consider the Higgsless case we set 
ch = in Eq. (A. 4) and consider the y/s <C mh limit using = \/2Xv. We find 

<Whi gg s oc \M\ 2 = ^0 + t) 2 (A.10) 
wheras for the case of a light Higgs with anomalous couplings we find from Eq. (A.5) 

"anom- higgs ^ l-^'l — 

under the assumption 



(Tanom-higgs OC \Mf = 4 ( s + *) ( A - U ) 



C -f»< (A.12) 
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This is how we arrived at Eq. (2.3). For rrih ~ 100 GeV this is true for scattering at the 
TeV scale as long as ch£, > 1/10. At lower values of the anomalous coupling the dominant 
effect comes from interference effects proportional to s + 1 instead of (s + 1) 2 . Thus smaller 
values of the anomalous couplings have a qualitatively different energy behavior that of 
larger values, making them especially difficult to resolve. 

B. Overview of Jet Algorithms 

While a comprehensive review of jet algorithms is beyond the scope of this work (see [62] 
for a recent review), here we will provide a short overview so the reader can quickly gain 
intuition into subjet techniques. 

Jet algorithms can roughly be divided into two categories: cone algorithms, which 
function as cookie-cutters to stamp out jets from calorimeter cells, and sequential recombi- 
nation algorithms, which build up a jet by merging four-momenta one by one in a prescribed 
order. Here we will focus on recombination algorithms. 

Each of these algorithms functions by defining a distance measure between every pair 
of four-momenta and for each four-momenta individually: 

d l3 =mm(PTlP 2 Tj) (ff) 2 ' d iB = P 2 T n j (B.l) 

for jets i and j. If the smallest distance measure at a given stage in clustering is between 
two four-momenta they are merged, otherwise the four-momenta with the smallest diB is 
declared a jet and removed from the queue. 



Algortithm 


n 


Approximate clustering order 


kx 


1 


soft— > hard 


C/A 





near— >far (in y-(f>) 


anti-fc^ 


-1 


hard^ soft 



Table 7: Parameterization and approximate behavior of sequential recombination jet algorithms. 

The different sequential recombination algorithms are distinguished by value of n ap- 
pearing in Eq. (B.l). These values determine the clustering order, whether one clusters 
beginning with hard four-momenta, soft four-momenta, or by angle (see Table 7 and Fig. 9). 
For the subjet analysis at hand, where reconstructing the softer subjet is essential, we there- 
fore use the k? algorithm which begins by clustering softer jets, preventing them from being 
cannibalized by the harder subjet. 
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